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Abstract 



Phase-contrast microscopy is used to monitor the shapes of micron-scale 
fluid-phase phospholipid-bilayer vesicles in aqueous solution. At fixed tem- 
perature, each vesicle undergoes thermal shape fluctuations. We are able 
experimentally to characterize the thermal shape ensemble by digitizing the 
vesicle outline in real time and storing the time-sequence of images. Analysis 
of this ensemble using the area-difference-elasticity (ADE) model of vesicle 
shapes allows us to associate (map) each time-sequence to a point in the zero- 
temperature (shape) phase diagram. Changing the laboratory temperature 
modifies the control parameters (area, volume, etc.) of each vesicle, so it 
sweeps out a trajectory across the theoretical phase diagram. It is a nontriv- 
ial test of the ADE model to check that these trajectories remain confined to 
regions of the phase diagram where the corresponding shapes are locally sta- 
ble. In particular, we study the thermal trajectories of three prolate vesicles 
which, upon heating, experienced a mechanical instability leading to budding. 
We verify that the position of the observed instability and the geometry of the 
budded shape are in reasonable accord with the theoretical predictions. The 
inability of previous experiments to detect the "hidden" control parameters 
(relaxed area difference and spontaneous curvature) make this the first direct 
quantitative confrontation between vesicle-shape theory and experiment. 
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I. INTRODUCTION 



Micron-scale fluid-phase lipid-bilayer vesicles have been observed in recent years under 
controlled laboratory conditions [p]-[TT[| to exhibit many amusing and diverse shapes. At 
the same time, there is now a one-parameter theory of vesicle shapes, the so-called area- 
difference-elasticity (ADE) model |1^-|15[ , which appears to be qualitatively consistent with 
available experimental observations. It would be nice, however, to have a vesicle with an 
accurately measured shape and known parameters, to plug these parameters into the theory, 
to predict a shape, and to compare with the measured one. Up to this time, this has not 
been possible, and, indeed, there have been few (if any) direct quantitative confrontations 
between theory and experiment. 

The reasons for this unsatisfactory state of affairs have their origins in both theory and 
experiment. On the theoretical side, the principal models which have been proposed to de- 
scribe vesicle shapes []T2,rB-TS] all have in common the same catalogue of stationary-energy 



shapes. Thus, simple observation of a vesicle whose shape can be found in the catalogue, 
while evidence for the general validity of bending-energy models, does not distinguish one 
variant from another. In order to test the model, it is necessary to study more indirect 
and/or delicate issues such as stability (absolute and relative) or shape-change systematics 
under variation of control parameters. This has not often been done [|I^ for fundamental 
experimental reasons: First, there are two important vesicle parameters which can be modi- 
fied systematically in the lab but are not subject to direct measurement. One of these is the 
spontaneous curvature, Co, which measures the preferred radius of curvature of the relaxed 
bilayer, based on the different lipid composition of the two constituent monolayer leaves 
and/or the different aqueous environments inside and outside the vesicle. This parameter 
is presumably the same for all vesicles in a single homogeneous suspension. The other is 
the relaxed area difference, AAq, between the two leaves, based on the different number 
of lipid molecules which they contain and the long relaxation time for lipid exchange 



between them pl|j2^ . This parameter will in general vary from one vesicle to another in the 
same suspension, based on the (unknown) manner in which vesicle closure occured during 
preparation and on any interleaf "flip-flop" or intercalation events which may have occured 
subsequently. In addition, the vesicle shape, itself, is a significantly ambiguous quantity. 
At laboratory temperatures, all nonspherical vesicles undergo significant and unavoidable 
thermal fiuctuations. Thus, at any nonzero temperature T, experiment must characterize 
a thermal shape ensemble. A single "snapshot," such as has often been exhibited in the 
previous literature, cannot do this. Finally, experiments have not in practice probed the 
full, three-dimensional vesicle shape but at best a two-dimensional cut through it at the 
focal plane of the observing apparatus. 

It is the aim of this paper to show how to deal with all these problems in a serious 
manner for the first time. Using video phase contrast microscopy, we recorded for each 
vesicle and at each temperature long time-sequences of two-dimensional vesicle contours. 
We parameterized these images in terms of a set of shape amplitudes. We used the shape- 
amplitude time- sequences to construct a thermal ensemble, from which we extracted a set 
of thermal expectation values. Using this data, we show below how to associate each vesicle 
with a particular point in the phase diagram. In principle, information is left over after 
the mapping, so that a nontrivial confrontation between theory and experiment is possible. 
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In practice, available resolution limits what we can do; nevertheless, nontrivial checks are 
possible. 

Overall, the results are encouraging. The vesicles that we have located in ("mapped 
into") the theoretical phase diagram do, generally, end up in regions where they are pre- 
dicted to be locally stable and to have low energies. Furthermore, observed thermal trajec- 
tories exhibit shape instability close to (if not always exactly at) positions predicted by the 
theory. Finally, after the instability the shape is in reasonable agreement with theoretical 
expectations. 

Gravitational effects play an important role experimentally. In order to record long time- 
sequences of vesicle shapes, it is convenient to adjust the density of the exterior solution so 
that the vesicles have a small negative buoyancy and collect at the floor of the experimental 
cell, where they remain within the focal plane of the microscope for long periods. In addition, 
gravity orients the long axis of vesicles of prolate shape so that it stays in or near the focal 
plane. These are practical issues. On the conceptual side, whenever it has non-neutral 
buoyancy, a vesicle is subject to gravitational shape deformations. The importance of these 
deformations has only recently been recognized and was not considered in the analysis 
of earlier experiments. In this paper, we first perform the full analysis without including 
gravity. Then, we devote a separate section to the consideration of gravitational corrections. 
The upshot is that gravitational effects can be significant; however, in the region of the 
phase diagram upon which we focus attention, there is no change in the qualitatively good 
agreement between theory and experiment. 

The layout of the paper is as follows: Section |I| introduces the theoretical background 
necessary to analyse the experiments. Section |T| describes the experimental procedures. 
Section |^ explains how the analysis of the experimental shape contours was carried out. 
Section |V] sets forth our results using a pure ADE mapping and ignoring gravitational 
effects. Finally, in Section 0, we explore the effects of gravity. Section [VII] provides a final 
assessment and summary. 



II. BACKGROUND 



A. The area-difference-elasticity model 



In order to have a language for discussing the experiments, it will be useful to present here 
a summary of some principal features of the ADE model. Additional material is available 
elsewhere [p!2|- [T5|JT8|j2^j25[| . At mesoscopic length scales, larger than molecular sizes but 



smaller than the persistence length, the shape 5* of a fiuid-bilayer vesicle is controlled by an 



energy functional W[S] consisting of two parts. The first, due to Helfrich ||2^, measures the 
overall bending energy and is scaled by the bending modulus k. The second requires a brief 
explanation: Assuming fixed bilayer separation, D, the actual area difference between the 
two leaves of the bilayer is 



AA[S] ='2D f dAH{v) 



(1) 



where H{v) is the local mean curvature at the point r of the vesicle surface and the integral 
runs over the (closed) vesicle surface. On the other hand, the preferred or relaxed area 
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difference, 



AAo= {Kut-N,^)ao{T), (2) 

of the two leaves is determined by the difference {Nout — Nin) between the number of hpid 
molecules in the outer and inner leaves. The relaxed area ao(T) per lipid molecule is a 
material parameter but can depend, of course, on the temperature T. Once the vesicle 
has closed, AAq can only change due to lipid fiip-fiop between the two leaves and/or lipid 
interchange with the aqueous environment of the vesicle, processes which are believed to be 
slow on the timescale of the mechanical shape changes we shall be discussing [0,^. The 
second contribution to W[S] measures the elastic energy necessary to force Ay4[S'] to differ 
from AAq, when the vesicle assumes the shape S. Because the vesicle is fluid, this (local) 
elastic strain is distributed uniformly over the vesicle surface and appears as an apparently 
nonlocal term controlled by a so-called nonlocal bending modulus R. The moduli k and k 
are both of order KD^, where K is the area stretching modulus of the bilayer so the 
ratio, 

a = k/k, (3) 

is generically of order unity. The material parameters k and R can be measured directly. 
For SOPC, it is believed that k ~ 0.90 ± 0.06 x 10"^^ J pi. It has been estimated that 



a ~ 1.4 |12|. (A somewhat higher value, k ~ 1.20 ± 0.17 x 10 ^^J and a comparable but 



quite uncertain value of a have been recently observed in tether-pulling experiments [^.) 
The energy scale k is much smaller than the energies necessary to change significantly the 
area A and volume V of the vesicle so these quantities may be regarded as fixed in 



comparing the energies W[S] of different shapes. 

Combining the two terms described in the previous paragraph (and dropping an irrele- 
vant, shape- independent term) leads to. 



G[S] + -{m[S]-mo 



(4) 



W[S] = 
where 

G[S] = ydA{2Hr, (5) 

which is the starting point of our theoretical discussion of shapes. In writing Eq. (^, we 
have chosen to rescale all lengths in terms of an "area length" Ra defined hj A = AttR\. 
Thus, the area difference appears in the reduced form, 

m[S] = AA/2DRa, (6) 

and the relaxed area difference combines with the spontaneous curvature Cq into a single 
effective reduced area difference, 

TfiQ = mo + 2co/a, (7) 

where mo = AAq/2DRa and cq = CqRa is the reduced value of the spontaneous curvature. 
Because Cq and AAq appear only in the combination mo, it is impossible in principle to 
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detect either one separately by a single shape measurement. Note that the bracketed terms 
in Eq. are all dimensionless ratios of lengths, invariant under a scale change of the shape 
S, provided that at the same time Cq is changed to keep cq fixed. In this sense, W[S] depends 
on the shape of S but not its overall size. An appropriate scale-independent volume measure 
is the reduced volume v = 3V/AttR\, which lies in the interval [0, 1]. 

To be a mechanically viable shape for a vesicle with given A and V, S must make the 
energy stationary at the corresponding values of v and rho [^, i.e., it must satisfy 

6W = = k[6G[S] - a{mo - m[S])6m[S]). (8) 

In general, there are several distinct branches of stationary shapes, which we label 
, mo), with corresponding energies {v,mo). To be a candidate for observation 
in the lab, a shape S^'^\v,Tho) must, in addition, be locally stable to small shape perturba- 



tions UsTJ. The lowest-energy branch (which must, of course, be stable) defines the ground 
state and should, in principle, be observed at sufficiently long times when the temperature 
is low. However, when energy barriers are large on the scale of ksT, other low-lying locally 
stable branches may remain metastable for long periods. For SOPC, k ~ 20 ksTroom p7[| , 
so metastability is expected to be common. 

Finally, we shall need below an important connection between the ADE-model shapes 



and those of the so-called spontaneous curvature (SC) model p6| , |32| , defined by the energy 
functional, 



SC 



[S] = '^j>dA{2H - Cof = k[G[S] - 2com[S] + const.), (9) 



which describes a model without differential area elasticity and having a spontaneous curva- 
ture Cq (cq is the corresponding reduced spontaneous curvature). The variation of Eq. (^ 
gives a condition which has the same form as Eq. (^ only with the replacement, 

2co = a(mo -m[^(")]). (10) 

It follows that any stationary shape S"*^"-* of the ADE model (^ is also a stationary shape of a 
spontaneous curvature model (^ with the cq defined by Eq. ([l^), which we shall henceforth 
refer to as the effective reduced spontaneous curvature for the ADE shape S*^"-*. Notice that 
the control parameters enter the variational shape equation (§) entirely via the coefficient 
of the second term. It follows that the stationary shape S^^^ is specified completely and in 
a way that is independent of a by giving v and Cq. For this reason, it will sometimes be 
convenient in what follows to think of the stationary shapes as S^'^\v,co) rather than as 
S^'^\v,fhQ), which still depends implicitly on a. The variation (P) may be thought of as 
proceeding in two steps: First make G[S] stationary at fixed m, thus defining a function 
(w, m), then subsequently carry out the variation with respect to m. It follows from Eq. 
(1) that Co can be evaluated as 

2co = = 2co + a(mo - m[S^-^]). (11) 

om 
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B. T = Phase diagram and the stability of prolate shapes 



The map in the {v, mo) plane of the regions where various branches provide the lowest 
energy shape constitutes the T = (shape) phase diagram for the ADE model. This phase 
diagram, which depends on a, can now be constructed rather easily numerically, at least, 
for V not too small and when mo/47r is not too far from unity P^J5B|J5^ . The experiments 
described in this paper deal with a branch of axisymmetric shapes, called "prolates," be- 
cause they have up/down symmetry and resemble prolate ellipses when they have reduced 
volume not too much below unity. The region of the ADE phase diagram in which these 
prolates appear is shown in Fig. |I|. It is bounded below by oblate axisymmetric shapes 
and by a region of non-axisymmetric shapes, which need not concern us further here . 



Above the prolates lies a region of pear shapes, for which axisymmetry remains but the 
up/down symmetry has been broken. This pear region is, in turn, bounded above by a line 
j^pear q£ fully "vesiculatcd" limiting shapes, consisting of two spheres joined by a narrow 
neck. The region above L^'^"-^ is incompletely explored but believed to be dominated by 
additional interconnected vesiculated shapes. The boundary between the prolates and pears 
at relatively high reduced volume involves a discontinuous shape change (corresponding to 
a simple crossing of the energy branches W^'^^^ and W^^") along the line JJp^o/p^"-"^ but a 
continuous shape change (corresponding to a bifurcation of W^^"""^ away from W^™) along 
the line labeled (7P''o/pear^ lower reduced volumes beyond the tricritical point T. Both 
transitions are often called "budding." 

It is important to emphasize that the prolate shape branch continues to exist outside of 
the "prolate region" of the phase diagram. Indeed, within the context of the ADE model 
and in the region of reduced volume shown in Fig. a stationary prolate shape exists 



for every value of mg These shapes can only be observed, of course, when they are 

locally stable. The region of local prolate stability includes the "prolate region" of the phase 
diagram but extends beyond it into metastable regions, where the true ground state has 
some other shape. It is a crucial test of the theory that prolate shapes observed in the lab 
should, indeed, map to the region of predicted prolate (local) stability. 

Metastability boundaries are marked by the first appearance of a soft mode, i.e., a family 
of fluctuations which lower the overall energy. The region of Fig. [1| within which prolates 
are predicted by the ADE model to be locally stable is bounded above by the line Mg™ and 
below by the line Mf +. These lines are calculated by an analysis of constrained Gaussian 



fluctuations about the calculated stationary shape |Q. The subscripts label the rotation 
mode |m| and (even/odd) parity of the sector where the first instability occurs. It is an 
important result of this theory that instabilities in sectors which break the symmetry 



are a property of the shape S''-"^ alone, while those in non-symmetry-breaking sectors depend 
in addition independently on a. For the prolate shapes, the boundaries C^'''°/p^"-^ and Mq™ in 
Fig. |l| both reflect instability in the symmetry-breaking sector (|m| =0, odd parity ). These 
boundaries are, thus, independent of a in a (f , cq) representation of the phase diagram, so, 
in the usual (f , mg) representation, they shift with a according to 

mo = m[^P™(t;,Co)] + ^. (12) 
These lines of shape instability (and not the actual shape (phase) boundaries!) are the 
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experimentally relevant (observable) ones. 

Strictly speaking, the above picture holds only in the low-temperature limit, since for 
T > sharp phase boundaries do not exist, because the vesicle is a finite system with 
finite energy and always explores its full phase space. Nevertheless, in practice, as long as 
the prolate branch remains locally stable and surrounded by energy barriers appreciably 
larger than ksT, thermally fiuctuating prolate shapes are readily seen in the lab. In the 
prolate region of the phase diagram, these fiuctuating shapes constitute a true, stationary, 
equilibrium ensemble. In regions which are only metastable, the set of prolate shapes should 
be regarded as a restricted ensemble, which may, however, be quasi-stationary for appreciable 
periods of time. In practice, metastability is expected to break down slightly inside the 
boundaries M^^°, when the metastability barrier becomes comparable to ksT. 



C. Prolate shapes for reduced volumes near unity: The hierarchy and the mapping 

We review briefiy here what is known theoretically about the T = stationary shapes 
S^^°{v, Co) for the relatively high reduced volumes which will be relevant for the experiments 
0,11. 

Prolate shapes are axisymmetric. Therefore, they are completely described by the curve 
made by their intersection with any plane which includes the symmetry axis. This curve 
may be written in terms of an arclength s which starts at the north pole {s = 0) and ends at 
the south pole (s = s*). We take the direction of the polar axis to be y and the perpendicular 
direction to be x. A representation which will be convenient for our purposes is 

^(^) = + E 4°^ sin(nvr-) , (13) 

^ n=l ^ 

where ip{s) is the angle between y and the outward-pointing normal to the curve. The 
overall length scale is set by s*. Note that ip{0) = and ijj{s*) = vr. The first term 
on the right describes a semicircular spherical vesicle shape. The coefficients 

{a^°^} parameterize deviations from the sphere. For shapes like prolates, which are up-down 
symmetric, the odd-ra coefficients vanish. We note as an aside that the coefficients {a^°^} 
cannot be set independently, since closure of the curve at the south pole requires that 

x{s*) = r ds cosi:{s) = . (14) 

<^ 

This places a complicated nonlinear condition on the set {ai"-*}, which for any real vesicle 
shape will automatically be satisfied. 

The stationary shapes of the prolate branch are given by the coefficients 02°^ (f, Cq). It is 
clear that, when v is near unity, the coefficients a^2n '^iH be small. It is a consequence of 
the stationarity condition (^ that these non-vanishing coefficients have the structure of a 
well-defined hierarchy ||25|| , 

A2{l-vf2 + B2{co) (l-v) + 0{{l-v)l) , 

B,{co) {l-v) + 0{ {l-v)l ) , (15) 
0( (l-^)i), 
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where A2 = (135/64)^/^ ^ 1.45 and the coefficients B2 and are hnear functions of Cq with 
coefficients of order unity. These results follow from Refs. |]T2| and [E^. It is a consequence 



of this structure that for v near unity 02°'' ^ ai"'* ^ ag°^ . . ., a hierarchy which reflects 
the fact that modes of higher n correspond to shorter wavelengths and cost more bending 
energy, is independent of cq at lowest order, because this is the only contribution of 
order (1 — f)2 and is, therefore, entirely determined by the constraint on the reduced volume 

V. 

Figure H, which was calculated by solving numerically the variational equation (H), illus- 
trates the dependence of af^ on v and cq for the prolate branch near v = 1. It is clear that 

3 

for V < 0.95 the terms of order (1 — v)^ and higher have an appreciable effect. Note that, 



except near the anomalous point |^ v ~ 0.85, knowledge of af' and v uniquely determines 



Co. 

This brings us, finally, to the issue of the "mapping" , i.e., of associating an experimentally 
observed T = prolate vesicle shape with a point in the ADE phase diagram. Fig. |l|. 
Knowing the vesicle shape means that we have direct experimental access to "geometrical" 
quantities such as v and m, through the shape coefficients a'^^. The abscissa, v, of the phase 
diagram is geometrical; however, the ordinate, rho, Eq. (0), encodes information about the 
initial area difference AAq and the spontaneous curvature Co, which are neither geometric 
nor directly observable in any other way. The solution to this apparent impasse is to use 
04 and Fig. |^ to infer a value of co{v,af^), the effective reduced spontaneous curvature, 
which is not observable, and to combine this with the then (theoretically) fixed m{v, co) to 
calculate mo (Eq. (|T0|)). Note that, in principle, any of the nonzero coefficients {af^} could 
be used to produce such a mapping, co{v,afn)- In practice, however, as the hierarchy (piSj) 
shows, a^2^ is very insensitive to co (because A2 is independent of cq and B2 is only weakly 
dependent on it) and Og^"* is sufficiently small so that experimental noise makes it a poor 
candidate. 

This framework is still incomplete in three senses. First, what the experiment observes 
is not a single T = shape but an ensemble of thermally fluctuating shapes. Second, the 
mapping as described above simply takes an experimental point and associates it in a one-to- 
one manner with a point in the theoretical phase diagram. It does not yet in any obvious way 
test the correspondence between theory and experiment. Third, effects of gravity should be 
taken into account. We discuss these important points in Sees. |V|, |V|, and [V^, respectively. 
But, before doing so, we turn to the experiments. 



III. EXPERIMENTAL TECHNIQUES 



A. Materials and Preparation 



For all experiments, vesicles were prepared from the common phospholipids 1- 
Stearoyl-2-01eoyl-sn-Glycero-3-Phosphatidylcholine (SOPC) or 1,2-Dimyristoyl-sn-Glycero- 
3-Phosphatidylcholine (DMPC). These lipids have their main phase transitions at 5° C and 
23° C, respectively [^]. They were purchased in powder form (Avanti Polar Lipids, Birm- 
ingham, AL, USA) and stored dissolved in chloroform: methanol (2:1) in special chemically 
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inert glass vials (Fisher Scientific) below -15° C. 

Preparation was done using standard techniques ||25|j40| . |4l[] : A few drops (30 fiV) of lipid 
solution (10 mg per ml chloroform: methanol) are spread with a syringe needle on a roughened 
Teflon disk. The solvent is evaporated in a vacuum chamber overnight. The disk with the 
dried lipid is placed in a glass beaker (50 ml) and pre-hydrated with a stream of Argon 
saturated with water vapor for about 20 minutes. Then, the desired solution for vesicle 
swelling is added and the beaker is covered with Parafilm and placed in the oven. To avoid 
heat shock, the solution and the beaker with the Teflon disk are heated separately to the 
swelling temperature prior to incubation. Swelling was done with 50 mMol sucrose solution 
at a temperature of 36° C. 

Successful vesicle development is indicated by whitish streaks in the swelling solution. 
These streaks are collected into an Eppendorf tube with cleaned glass pipettes and incubated 
at the swelling temperature. Excess glucose solution (48 mMol) is then added to obtain the 
desired density for the vesicles in the observation chamber. The end result of this procedure 
is a vesicle suspension with an interior sucrose solution and an exterior glucose solution (with 
a slight admixture of sucrose). The excess density of the interior relative to the exterior sugar 
solution is approximately 3.3 g/1. This is needed in order that the vesicles sink gently to 
the bottom of the experimental cell, as discussed further below. Vesicles were stored at the 
swelling temperature and used within a few days. 



B. Experimental Setup and Data Acquisition 

For observation, vesicles are placed in a specially designed microchamber, tightly sealed 



with glass plates above and below to prevent evaporation ||2^. Temperature is monitored 
by a thermocouple inserted into the observation chamber. A water bath, incorporated 
integrally into the chamber, provides temperature uniformity and control at the level of 
±0.1° C. Because their density is slightly higher than that of the surrounding solution, the 
suspended vesicles fall to the bottom of the cell, where they rest gently against the lower 
plate and are observed from below via video phase-contrast microscopy. We use a standard 
inverted Leitz microscope equipped with phase contrast (Leitz Phaco 40/0.65), capable of an 
overall magnification of 500 times and incorporating annular illumination and a phase ring. 
The light source was a Hg arclamp powered by a high-voltage transformer. A permanent 
green filter and various grey filters were used to minimize degradation of the lipids. The 
video camera was positioned above the eye piece in such a way as to gain a resolution of 86 
nm per pixel in a 480x480 frame. 



The visible phospholipid structures which collect at the bottom of the observation 
chamber in the microscope's focal plane are typically very diverse including topo- 

logically complex and multilamellar structures, small vesicles included within larger ones, 
vesicles connected to one another by sub-microscopic tethers or tubes, vesicles with obvious 
adhesions, etc. For detailed observation, we try to select simple, topologically spherical, unil- 
amellar structures, without identifiable microscopic connections or adhesions. In addition, 
we monitor the fluctuations of each candidate vesicle for some time prior to data acquisition 
in order to reject those with obviously "abnormal" behavior, e.g., those exhibiting unex- 
plainable asymmetry or sudden changes in apparent area or volume. The final fraction of 
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usable vesicles is less than one percent. 

Images of selected vesicles are simultaneously displayed on the video monitor, saved to 
tape (U-Matic, Sony), and processed in real time, as will be described in the next subsection. 
To sample at a fixed temperature a single thermal "shape ensemble" takes about 20 min. 
Recording a thermal shape trajectory requires data at several different temperatures for the 
same vesicle. The chamber is allowed to equilibrate for at least 5 min. after each temperature 
change. Temperature is recorded with a precision of ±0.1° C. The total amount of data 
gathered consists of over 80 hours of video tape of more than 150 vesicles, including a wide 
range of shapes. Budding is an ubiquitous process, which we observed at least 15 times in 
a controlled fashion. It is important to get long runs at each fixed temperature in order 
properly to sample the full thermal shape ensemble. On the other hand, thermal-trajectory 
runs which extend over more than a few hours appear to be contaminated by systematic 
drifts, presumably due to lipid degradation and/or fiip-fiop between bilayer leaves, which 
establishes intrinsic limits on sampling density and run time. 

In this paper we restrict analysis to three particular SOPC vesicles (A, B, and C), 
which have in common that they started with prolate elliptical shapes and, on heating, 
eventually underwent budding transitions, as illustrated in Figs. 0-^. Qualitatively, other 
vesicles monitored behaved similarly, although they followed different trajectories, some 
exhibiting sharp shape transitions and others not. The reason for selecting the budding 
trajectories is that the location of the budding instability provides a particularly stringent 
test of the theory, as we shall discuss in Sec. 0. 



C. Processing the Video Image 

In order to analyse the data, it is necessary to reduce the video image to a time sequence of 
digitized shape contours. This was done by using a frame grabber (Matrox, Dorval, Quebec, 
Canada) to capture each image, computer processing the image in real time, storing the 
digitized contour point in memory, and then grabbing a new image. The image processing 
algorithm (described below) requires between 0.4 and 0.6 s (depending on vesicle size) on 
a PC with a 486 DX CPU and a 66 MHz clock speed. Thus, for a video frequency of 30 
frames/s, we are processing every 15th frame. This is relatively slow compared to processing 
times on the order of 0.1 s, which have been reported in the literature P,P3|. However, in 
contrast to these fast procedures, our algorithm has a better-than-pixel accuracy in finding 
the contour P^J15[ |. This high resolution turns out to be critical to the success of our 



experiments, since we shall need to resolve small changes in mean shape in tracking the 
thermal trajectory HH . 



In phase-contrast microscopy the image of the vesicle edge exhibits a "halo," with a light 
band (intensity maximum) just outside the vesicle and a dark band (intensity minimum) 
just inside. Typically, the intensity profile crosses the gray value of the local background at 
its steepest point, and we have taken this point to be the nominal position of the vesicle 
boundary P7| . 

The contour-digitizing algorithm is fully described in Ref. p5|. The algorithm works 



on the (integer) pixel grid {nx,ny) and requires initialization by hand to the vicinity of the 
outline of the particular vesicle to be studied (there are ordinarily several vesicles in the 
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field of view). Suppose that a scan in at fixed = n'^'^ intersects the halo. By averaging 
the grey values in the vicinity of the halo, we establish a local background intensity. The 
profile of intensity-versus-n^. crosses this background value at a point x (generally non- 
integer) which may be determined by linear interpolation. The point (x, n^°^) is then stored 
as a contour point, and the algorithm steps Uy ^ Uy + 1 and starts again. Note that the 
interpolation procedure allows determination of the x coordinate of the contour point with 
better-than-pixel precision. Whenever the contour profile becomes steeper when scanned 
in the y direction rather than in the x direction, the algorithm automatically switches to 
scanning at fixed n^., and vice versa. Each contour is terminated at closure. Motion 
of the vesicle between successive frames is normally small enough so that each subsequent 
frame can be started where the earlier one terminated, so the initialization step needs to be 
carried out only at the beginning of each run 

The digitized contours exhibit noise at the pixel level (1 pixel=86 nm). This behavior 
presumably reflects the intrinsic noise of the original optical signal, the pixelation statistics, 
the digitization of the grey scale, and other factors. To remove some of this microscopic noise 
before data analysis, it is convenient to smooth the observed contours. This was done by 
applying a tenth-order binomial fllter to the x and y contour coordinates, thus averaging 
over an effective width of about 5 pixels. The distribution of deviations of the original data 
points from the smoothed contour is Gaussian with a typical full width at half maximum of 
about 0.7 pixels, thus giving an effective local lateral resolution of about 30 nm (compared 
to a typical vesicle size of several microns). This resolution, well below the nominal optical 
resolution given by the wavelength of light, illustrates the delicate line-shape discrimination 
achievable via phase contrast ||5^ and is more than adequate for quantifying the overall 
vesicle shape and the low-lying fluctuation modes. 

The result of this process is a time-sequence of several thousand digitized contours, 
illustrating the shape ensemble of each vesicle at each temperature. The relationship of these 
two-dimensional contours to the three-dimensional vesicle shape requires a brief discussion. 
The "general wisdom" seems to be that what is seen in phase contrast microscopy is a cut 
through the vesicle in the focal plane P, |i3| , |51| -p^ . However, this is an over simpliflcat ion. 
Phase contrast is particularly sensitive to edges, so vesicle boundaries which "overhang" the 
focal plane (relative to the optical axis) may contribute to the image to a greater or lesser 
extent depending on the focal depth and the amplitude of the edge contrast. Following the 
practice of the literature, we shall ignore such effects in what follows. We wish only to point 
out that there are substantive issues here which deserve to be addressed more fully in future 
work. 

In collecting data, the microscope is focussed on the maximal cross section of the vesicle 
under observation, and this focal plane does not change over time. For vesicles such as 
A, B, and C which are (on the average) prolate and axisymmetric, gravity tends to orient 
the symmetry axis horizontally, i.e., to bring it into the focal plane. Thus, the contours 
(such as those based on Figs. 0-§|), which constitute our raw data, may be thought of as an 
ensemble of cuts through the mean symmetry axis of the fluctuating vesicle. Fluctuations of 
the symmetry axis out of the horizontal plane modify this simple picture: If the focal plane 
no longer includes the symmetry axis, then the depth of focus and the edge enhancement 
mentioned in the previous paragraph probably give the resulting image the character of a 
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projected outline of the tilted vesicle. In principle, this projected shape is different from a 
true axial section. In practice, the stabilizing effect of gravity is large enough so that these 
out-of-plane fluctuation effects are almost always small (except near the spinodal line of 
the budding transition), so we will treat the two-dimensional contours as if they represent 
axial sections. Note that there is a balance here. In order to keep the theoretical analysis 
simple, we would like to ignore the effects of gravity on the vesicle shape. On the other 
hand, in order to perform the experiment conveniently, we use gravity to localize the vesicle 
in the bottom of the chamber and to orient the symmetry axis (of prolate vesicles) to the 
horizontal plane (see further discussion at the end of Sec. |Vp. 

IV. ANALYSIS OF THE TWO-DIMENSIONAL SHAPE CONTOURS 

This Section describes how we parameterize the individual digitized two-dimensional 
shape contours discussed in Sec. how we average the shape parameters over each thermal 
ensemble, and how we infer T = three-dimensional shape information from these averaged 
parameters. 



A. Thermal fluctuations of the vesicle shape: General discussion 

Varying the temperature has two effects. On the one hand, it modifies, through ordi- 
nary thermal expansion, the temperature-dependent control parameters, A{T), V{T), and 
AAq{T), as well as the material control parameters, k, R, and Co, which appear in the 
Hamiltonian (p. These effects produce the so-called "thermal trajectories," which we shall 
discuss in Sec. |V|. On the other hand, even if all these parameters were temperature in- 
dependent, there would still be ordinary thermal fluctuations. It is for the moment these 
purely thermal fluctuations to which we direct our attention. 

There is, in principle, no way of taking a single fluctuating shape contour and inferring 
the corresponding T = shape. At best, we must take a full thermal shape ensemble and 
use theory to infer the T = shape of the vesicle with the same control parameters. When 
the fluctuations are large, even this is beyond present theoretical capability. However, when 
fluctuations are small enough so that they may be treated at the Gaussian level, progress 
can be made. 

The upshot of a recent study of Gaussian fluctuations of vesicles of arbitrary axisym- 
metric shapes |]37| , |3^ , |55| may be summarized as follows: Any typical fluctuating shape may 



be regarded as a T = shape appropriately translated and rotated (the so-called Euclidean 
modes) plus an area- and volume- conserving normal (i.e., perpendicular) displacement u{r) 
at each point r of the surface. It is a special feature of these fluctuations that (because 
of the strict area and volume constraints) both the average displacement {u{r)) and the 
mean-square fluctuations {u^{r)) are generically of the order ksT/u, where u is a. typical 
static fluctuation-mode energy. Note that the rms fluctuations are always larger than the 
shift when the fluctuations are small. In the analysis which follows, we shall assume that the 
Gaussian regime holds and we shall ignore the mean shifts. Ordinarily, the fluctuation-mode 
energies u are of the order k. For our vesicles K/ksT ~ 20, so for most regions of the phase 
diagram, this is an excellent approximation. There is, however, an important exception. At 
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the instability boundaries (Fig. |l|) one of the modes becomes soft. Thus, near enough such 
boundaries, the Gaussian treatment fails, and we may expect difficulties (see Sec. 0). 



B. Parameterization of the two-dimensional shape contours 

We interpret the measured two-dimensional contours as being sections which include the 
principal symmetry axis of the vesicle (Sec. [Ill C| ). Thus, in the spirit of the last paragraph 
and up to corrections which are normally of order ksT / the center of mass is located 
in (i.e., near) the focal plane at (i.e., near) the point which is the center of mass of the 
digitized contour. We determine this point numerically for each contour. We then find the 
(approximate) principal (long) vesicle axis by calculating the two-dimensional moment-of- 
inertia tensor with respect to the center of mass and diagonalizing. This determines the 
(nominal) principal axis of each vesicle shape and identifies the north and south poles. We 
call the direction of the principal axis y and the corresponding perpendicular direction x 
(which is, of course, not necessarily a principal axis of the three dimensional vesicle). Thus, 
each experimental contour is reduced to a set of points In what follows, we treat 

each half-contour separately (each image has two half contours) and take Xj > 0. 

It is convenient to represent each half-contour in the angle-arclength ip{s) representation 
of Sec. [11 C| by calculating 



ipi = — arctan(^^^^!^^^ — ^^-J-^ ^ (16) 

where the arctangent is defined on its Riemann surface, i.e., ip{s) is continuous at the equator. 
The arclength s is measured from the north pole. The parameterization parallels Eq. (|T^, 



s 



V^(-s) = vr— + ^ a„sin(n7r— ) . (17) 



s ^1 s 



The coefficients {a„} are obtained by a numerical integration using the trapezoidal rule, 

an = {-IT- + J2i sin(n7r^) + ^i+i sin(n7r^) ) (18) 
n ~^ s s s 

and, henceforth, they replace the points {xi,yi} in representing the half-contour. (M is 
the number of digitized points in the half-contour.) Note that the contours here are not 
up-down symmetric, so the odd-n coefficients do not in general vanish, as they did for the 
T = prolate shapes. Similarly, ^/'(O) and ipi^*) are normally nonzero. 

C. Thermal ensembles and T = shapes 

For each half-contour of each video image, we calculate the shape coefficients {a„} plus 
the nominal ( "effective" ) vesicle area and volume, 

M 

Ae = irY^ixi + Xi+i){si+i - Si) (19) 

i=l 
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and 



= -X E(^' + ^li)iyi+i - Vi)- (20) 



In a similar spirit, we compute an effective reduced volume, 

Ve 



(21) 



For axisymmetric vesicle shapes, these equations would calculate the true area, volume, and 
reduced volume, respectively. Since each image is only a snapshot of a section through a 
fluctuating vesicle, is only approximately equal to the true reduced volume v. Note that 
Ve fluctuates in time for successive images of a given vesicle at a fixed temperature, while 
V is in principle constant, since the true area and volume are conserved during the shape 
fluctuation. 

The several thousand images which constitute a typical experimental run with a given 
vesicle at a fixed temperature lead to characteristic time series for the quantities and 
Ue, as illustrated by Figs. |^ and ||. Although the series are noisy, we expect to see memory 
effects between successive images, as long as there are any characteristic physical relaxation 
times longer than the 0.5 s between successive grabbed images. We have estimated else- 
where |TD| the typical relaxation times expected for these vesicles. Away from instabilities, 
the characteristic times are expected to be at most several seconds, which is consistent with 
direct visual observations of the optical image. As the vesicle approaches an instability 
(which occured at f = 0.878 for vesicle A), one sees rapidly increasing relaxation times, cor- 
responding to a spinodal slowing-down JlOl . This tendency is clearly visible when comparing 
Figs. 1^ and ^ As long as the data set spans a time interval much larger than the longest 
relaxation time, we may expect that the time sequence samples an effective stationary en- 



semble [Q. In this sense, the thermal ensemble of fluctuating shapes is characterized by 
the set of ensemble averages (a„), (ama„), (fe), (fg), etc. Indeed, the distribution functions 
P(a„), -P(fe), etc., are typically Gaussian form in shape PB[ . 



In order to proceed with the mapping, we now need a procedure for inferring the zero- 
temperature quantities v and {a^°-*} from the thermal data. We do this in the crudest way, 
by simply making the identifications, 

V = {v,) and = (a„) . (22) 

As explained above (Sec.[iV A]), the justification for these identifications is that the averaging 



process suppresses the rms fluctuations, which are of order y /c^T/u;. This leaves the thermal 
shifts plus the terms of order (u^), both of which scale as ksT /u, which we ignore in first 
approximation These values of v and 04°'' allow us (Sec. |II C|) to infer cq (and, thereby, 
mo) from Fig. ^, and, thus, to complete the mapping. 

It is hard to give any meaningful estimate of the real uncertainty in the derived quantities 
V and a^^\ For a truly stationary ensemble, the purely statistical (sampling) uncertainties 
in the average quantities should decrease as the ensemble sampling becomes denser. In 
practice, our runs are necessarily of finite length (Sec. [Ill Bp . Indeed, if we divide the data 
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set into two parts, corresponding to earlier and later times, we typically see a spread of 
values corresponding typically to ± 0.001 for v and ± 0.002 for the {a^^-'j's (and somewhat 
larger near the budding instability). It is this measure which we adopt as an estimate of 
the statistical uncertainties. Of course, there are also systematic errors, such as the thermal 
shifts (which we have neglected) , the failure of the Gaussian picture (where fluctuations are 
large), the sampling error (when the relaxation times are long), the fluctuations of the major 
axis out of the focal plane, and the effects of gravity (Sec. |V^. The upshot is that well away 
from the instability boundary the statistical uncertainties are probably realistic, except for 
the systematic influence of gravitational effects. Near the instabilities, the situation is less 
well defined. These statistical uncertainties translate (via Fig. |^) into uncertainties in cq, as 
we shall illustrate in Sec. 0. 



V. RESULTS (WITHOUT GRAVITATIONAL CORRECTIONS) 

Each of the three budding vesicles, A,B, and C, started at a relatively low temperature 
with a nearly spherical shape (i.e., v ~ 1). As the temperature was raised, the reduced 
volume decreased, until at a certain temperature (different for the different vesicles) a "bud- 
ding" instability occurred (see Figs. ^ and|^), i.e., the vesicle suddenly necked down and, 
over a time interval of 1 — 10 seconds, developed a small quasispherical satellite. (This 
time range is due to the different vesicle sizes, since typical relaxation times scale with the 
third power of the vesicle radius [0.) Up to the budding threshold, the thermally induced 
changes in the fluctuating ensemble are reversible to within experimental precision. The 
budding process, itself, is a mechanical instability [|I^. In fact, the budding can be reversed, 
but only by cooling to a temperature significantly below the budding temperature |^ . The 



size of the fluctuations and the scale of the longest relaxation time increases dramatically 
as the temperature approaches the budding temperature (see Figs. ^, H, 0, and H). These 
effects have been interpreted in terms of a simple Landau theory [0. 

Our results for the three budding vesicles. A, B, and C, are summarized in Table |I[ The 
average amplitudes (a„) were generally very small for odd n, as expected in the prolate 
phase [^]. The even coefficients (ag) and above were too small to distinguish from zero, 
presumably because of the hierarchy (|15]). The value Vb of the reduced volume at budding 
was determined by extrapolating the experimental temperature dependence v{T) to the 
observed budding temperature Tf,. 

Figure |^ shows the result of mapping this data into the theoretical (v, Cq) diagram by 
using (ve), (04), and Fig. ^, as explained in Sec. [II C| . The instability lines M'^'^"'^ and M°^'' 
are just the appropriately mapped versions of the corresponding spinodal lines of Fig. 
The advantage of this representation is that it is completely independent of the value of 
a, as explained after Eq. (0). Theory predicts that the prolate shapes are locally stable 
only between the two spinodals. With the exception of the highest-temperature point in the 
trajectory of vesicle A, we see that the mapped shapes do lie in this region. Fig. |l^ shows 
the same data plotted in the {v, mo) phase diagram. Fig. ^ The required relation between 
Co and fho is based on Eq. (0). In order to evaluate mlS^""^], we solve the variational shape 
equations derived from energy functional Eq. (01) for the given values of v and cq. We have 



taken a = 1.4 in making this transformation |28,32|. Since the spinodal boundaries map 
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right along with the data points, there is no change in the predicted stabihty. 

The fact that, with a single exception, the mapped data points lie neatly sandwiched 
in the region of predicted (local) stability is a stringent quantitative test of the theory and 
constitutes the single most important result of this work. We emphasize again (see Sec. 
p.1 B|) that the shapes mapped by Fig. ^ are variationally stationary (by construction) but 
not necessarily locally stable, so that an arbitrary shape could end up anywhere in the phase 
diagram. 

It is worth pointing out that fact that the values of Cq derived from the data points are 
all of order unity (as expected on the basis of the theory) is also an important test. Fig. |^ 
shows that Cq values between —5 and 10 are associated with values of of'' in the narrow 
range between —0.02 and 0.02. If the theory were significantly in error, it would be quite 
easy to have produced very large or very small values of cq . 

Indeed, in a certain sense, our ability reliably to distinguish shape changes corresponding 
to differences of order unity in cq is, in itself, surprising. Consider that, for a vesicle of radius 

10 iim. at a reduced volume v = 0.9, a difference in cq of ibl corresponds to a change in shape 
which modifies the pole-to-pole contour length s* by only 20 nm. This number (the smallness 
of which is a direct consequence of the hierarchy) is below the local lateral resolution of the 
contour. How is this possible? First, one has to realize that one does not measure a single 
distance only. Rather, the amplitudes are calculated globally from an integral (see Eq. ([T^)) 
over about 600 contour points, each of which deviates from the reference shape. Second, 
one is interested in a low mode, which is insensitive to local perturbations in the membrane. 
And, third, the amplitudes are averaged over typically several thousand contours, giving an 
effective sample size on the order of 10^. Thus, shape differences on the 10 nm scale are 
detectable 0. 

The "thermal trajectories" corresponding to each vesicle encode the effect of the experi- 
mental control parameter (temperature) on the quantities v and Cq (or mo); defined in Sec. 

11 A| . These quantities, in turn, depend on the volume V, area A, relaxed area difference 
AAq of the vesicle, on the thickness D and spontaneous curvature Cq of the membrane, 
and on the ratio a of elastic constants. All these quantities are in principle temperature 
dependent, and, if these dependences were known, we could calculate the thermal trajectory 
and compare with that found in Figs. |^ and |T^. The volume thermal expansion coefficient 
{Pv ~ 3 X 10~^/K for water) is known to be small compared to the area thermal expansion 
coefficient {(3a ~ 3 x 10~^/K for SOPC Q). It is also known that the total bilayer 
volume AD is only weakly temperature dependent. A simple model is to assume that only A 
and D are temperature dependent. When Cq = 0, as is reasonable for a symmetric bilayer, 
this assumption leads to the simple result that the product vuiq is temperature independent 
|T^,|23|. This hypothesis predicts trajectories of the general shape and scale shown in Figs. |^ 
and IC but significantly less steep than those observed. It is not hard to make more refined 
models consistent with the data, for example, by using a non-zero spontaneous curvature 
Co and/or a differential thermal expansion for the two leaves of the bilayer [p!6| , p5| . Unfor- 



tunately, direct measurements of these quantities are not available, so no useful conclusions 
can be drawn at this stage. 

Another set of evidence bearing on the consistency of the observations with the theory 
is the relative size and shape of the main vesicle and the bud which forms at the instability. 
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At the crudest level, theory predicts that the final state after budding will be pear shaped 
(rather than fully vesiculated, with a microscopically narrow neck) when the budding occurs 
for reduced volumes less than Vc = 0.875, where the spinodal crosses the limiting-pear line 
j^pear ^ggg pjg_ As tt increases, this crossing point moves to higher values of v, so the 
observation that vesicle A buds to a fully vesiculated state would be inconsistent with a value 
of a larger than 1.4. This observation places an upper bound on a. Another qualitative 
prediction is that the ratio r of the bud radius to the radius of the remaining main vesicle 
should increase as the reduced volume at budding decreases. Thus, r should be largest for 
vesicle A and smallest for vesicle C, as is, indeed, observed. On the other hand, on the 
limiting line L^^"-"^ ^ the vesiculated configuration consists of two spheres, so there is a simple 
relation between the ratio r and reduced volume v ||3^ , |62|| . In particular, close examination 
of Fig. I provides a value of r close to 2.8 (pure geometry, since the shape is very close 
to being two spheres). This corresponds to v{r) = 0.873. Indeed, budding is observed at 
Vb = 0.878 ± 0.002 for vesicle A. 

Overall, the agreement of theory and experiment is reasonable with the exception of the 
last point of the thermal trajectory of vesicle A, which lies distinctly above the theoretical 
spinodal boundary M^'"" (Figs. |^ and |1^). This last point is worrying. Indeed, even below the 
spinodal line, there should be a (fuzzy) unstable region where the energy barrier out of the 
metastable state is of order fc^T. (The fact that vesicles B and C appear to bud increasingly 
below M^*"" for higher reduced volume suggests that there may be some systematic effect at 
work which is distorting the locus of instability.) We have considered three possible reasons 
for this discrepancy. 

First, our identification (|22| ) involves the assumptions that (a) the fiuctuations are small 
enough to be treated at the Gaussian level and (b) the (Gaussian) thermal shifts and rms 
fiuctuations (of order ksT/u) can be neglected. At the spinodal, fiuctuations diverge [ITU 



so neither of these assumptions is valid, and the identification ( p2D is expected to fail. It is 
entirely plausible that these assumptions are already breaking down near the spinodal, at 
the last stable point. Because the effects of fluctuations beyond the Gaussian level have yet 
to be calculated, we cannot at this stage assess the impact that such corrections might have 
on the near-spinodal points of trajectory A 

Second, we have assumed that the major prolate axis is (effectively) in the focal plane of 
the microscope. If this axis is appreciably out of the focal plane, then the digitized images 
cannot be thought of as sampling axial sections of the three-dimensional fluctuating shape, 
and the whole analysis of these images would have to be redone. As long as fluctuations 
are small, it is reasonable to assume that gravity acts to keep the prolate axis aligned. But, 
near the spinodal line, there are large, slow pear-like fluctuations |Ty], which are not "up- 
down" symmetric (i.e., which break the symmetry between the north and south poles). In 
this situation, gravity may be expected to systematically re-orient the small end of the pear 
towards the bottom of the chamber, thus tipping the effective symmetry axis away from 
the horizontal. Once tipped, the symmetry axis is inhibited by gravity from returning to 
the horizontal, so one expects long intervals of asymmetric, pear-like data to appear in the 
near-spinodal time-sequences. In fact, the data for the last point of the vesicle A trajectory 
do show an "anomalous" cluster of frames with simultaneously large 03 and 04, and these 
frames exhibit a fuzzy contour proflle near the small end of the pear, indicating an overhang 
of the vesicle membrane beyond the focal plane. The effect of excluding this segment of the 
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time-sequence is to lower (04) nearly to the spinodal, thereby improving agreement between 
theory and experiment. 

Finally, we have so far treated the effect of gravity as only something which positions the 
vesicles at the bottom of the chamber and aligns (prolate) axes in the horizontal plane. In 
fact, it will also modify the zero-gravity shape analysis which has been up to this point the 
basis of our mapping procedure. What effect do gravitational shape changes have on the 
analysis and can they explain the observed discrepancies? This is the subject of the next 
Section. 



VI. GRAVITATIONAL EFFECTS 

A. Qualitative considerations 

When the density of the solution which fills the vesicle interior is greater than that of 
the exterior solvent, the vesicle will fall to the bottom of the container and, once in contact 
with the bottom, will deform in such a way as to decrease the gravitational potential energy 
of the interior, higher- density material. The overall shape involves a balance between the 
previous bending energy (0) and a new gravitational energy. 



Wgrav[S]=9oApJzdV, (23) 

where go is the local acceleration due to gravity, Ap is the excess mass density of the interior 
solution, z measures height above the bottom of the chamber, and the integral is over the 
interior volume of the vesicle. 

The ratio of the energy scale goApR\ of this gravitational term to the scale k of the 
bending energy defines a dimensionless parameter, 

J-^^^, (24) 



which measures the relative size of gravitational and bending energies p3[- When g is very 
small, we may expect shapes which are not significantly deformed relative to the gravity-free 
case. When g is very large, gravitational energy dominates and vesicles will tend towards 
circular pancakes [Q, squashed against the chamber bottom, insofar as constraints on area 
and volume allow. (Of course, if f = 1, then the vesicle can only be spherical.) 

Experimental values for our vesicles A, B, and C were nominally g =2.2, 0.3, and 1.9, 
respectively 0]. Here, we use a value of k = 0.9 x 10"^^ J for our estimation 0. We may, 
thus, expect gravitational corrections to be appreciable for vesicles A and C but relatively less 
important for vesicle B. The qualitative effect of gravity on the mapping is not hard to see. 
Roughly speaking, a "pancake" deformation will make the focal-plane section of a prolate 
rounder and larger in area than it would otherwise be. Therefore, gravitational corrections 
will lead to larger values of Vg. Although the coefficients a„ must approach zero for large g, 
it is not obvious where the asymptotic regime sets in, so the sign of the gravitational shift 
in a4 cannot be inferred a priori. In order to estimate these effects quantitatively, we need 
to be able to calculate vesicle shapes in the presence of gravity. 



18 



B. T = shapes in the presence of gravity 



The only previous calculation of vesicle shapes including gravitational effects was done by 
Kraus et al. [^. These authors found a gravity-induced prolate-oblate transition for values 
of g similar to those encountered in our experiments. Following this work, we use a polyhe- 
dral discretization of the vesicle surface and employ the program Surface Evolver to 
search iteratively for the shape-energy minimum. Numerical minimization in the presence of 
a hard-wall constraint for the chamber floor leads to special problems in stability. For this 
reason, we replaced the hard wall by a soft substrate potential, V„{z) = Vo^exp{—z/zo), 
with Vow = and zq = O.IRa- These parameters seem to provide a good compromise 
between numerical stability (favored by a softer potential) and a deformation of the shape 



caused by the soft tail of the potential which is as small as possible |^3[. Constraints on 
area and volume are respected. The energy is minimized by moving the vertices in the 
direction of the energy gradient or, alternatively, by a conjugate gradient method. Symme- 
tries such as mirror planes can be exploited. For most of this work, only a vertically-cut 
quarter section of the vesicle was actually computed. It turns out that the results for the 
final shape and energy are very sensitively dependent on the triangulation in a way that 
we cannot completely control. We have tried to overcome this problem by fitting a linear 
interpolation to a grid of data points, as described in Sec. VI B below. This procedure 
averages out random fluctuations from one point to another but cannot address any subtle 
systematic dependence on grid size which might be hidden beneath the fluctuations (we did 
check directly for such a grid-size effect, and none is apparent at the level of accuracy we can 
achieve). In the absence of a more reliable measure, we have simply used the deviation of 
the computed data points from the smoothed interpolation to give an estimate of the error 
introduced by the triangulation. 

Finally, we point out that, not only does gravity influence the shape and energy of a vesi- 
cle at given v and Cq, but it also changes the relative energy of different shape branches, thus 
shifting phase boundaries and stability boundaries in the phase diagram. Thus, in looking 
for gravitational corrections to the experiments, we also need to compute the gravitationally 
shifted spinodal line M^^". We have done this in the (w, cq) representation (Sec. [II A|) for 
g = 2.2, in which now (c.f., Eq.(|Tl|)) 

dG{v,m) ldWgrav{v,m) 

— dm K dm = ^Cq + a{mo - m[S^ >\). (25) 

This involves two numerical fits to Surface Evolver data, first for the computation of 
m{M^^'°) and then for the derivative which evaluates cq. As a consequence, the quality of 
the results is rather poor. As shown in Fig. |^, the result for our experiments is a shift 
of the upper spinodal upward in cq by about one unit; but, the numerical uncertainties 
are unfortunately comparable in size to the shift. More detailed calculations would require 
finer triangulations and much longer relaxation times. Since the computational investment 
is already appreciable and experimental uncertainties are already large near the spinodal, 
additional investment at this time does not seem wise. We did not compute the location of 
the full lower spinodal including gravity; however, for Cq = (and g = 2.2), we do know ||23 



that the prolate-oblate transition occurs at w = 0.94 (as plotted), which also corresponds to 
a small upward shift. 
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C. Gravitational corrections to the mapping 



We estimated gravitational effects by running the Surface Evolver program for v be- 
tween 0.875 and 0.975 in steps of 0.025 and for cq between 2.5 and 10.0 in steps of 2.5. For 
each pair of these parameters, we computed shapes over a range of small g values. For each 
shape, we took a maximal horizontal section and computed effective values of Ve-, a2, and 
04 using formulas (|T8|) and (|2l|) . These values varied in a roughly linear way with g, only 
with some superimposed fluctuations which we attributed to the triangulation sensitivity 
mentioned in the previous subsection. We then assigned effective values of Ve{g), 0.2(9), and 
a4{g) for the gravitationally distorted shapes by making a straight-line fit to these computed 
points passing through the values previously computed for g = 0. For the values of g, v, and 
Co relevant to the experiments, the gravitational shift of 04 is comparable to the numerical 
errors. On the other hand, the difference between Ve{g) and the actual reduced volume 
V (see Fig. |TTp does lead to a noticeable correction in the values of Cq inferred from the 
experimental data for vesicles A and C. 

This family of lines may then be used to calculate gravitational corrections to the ex- 
perimental data. Since g is known, it is only necessary to take the measured values of Vg 
and a4, which belong (presumably) to shapes which are gravitationally distorted, and to 
infer the corresponding values of v and Cq. Figure |n| shows, for example, the calculated 
gravitational shift in the apparent volume, (fg — v), at g = 2.2 (appropriate for vesicle A) 
for representative values of v and Cq. Note that the Ve is always larger than v, in agreement 
with the qualitative argument of Sec. VI A, so that gravitational correction always shifts 
the data points to the left in Fig. The corrections increase for values of cq close to the 
prolate-oblate transition, i.e., for small Cq, where increasing the volume of a prolate vesicle 
at fixed Co leads eventually to a transition to an oblate shape with a vertical symmetry axis, 
thus producing a circular focal-plane section (i.e., = 1). The behavior of the shift of the 
apparent volume near this prolate-oblate transition can be understood as follows: Above 
the transition, for v > Vc{co,g), the focal cut of the oblate vesicle is circular, independent of 
(co,^), i.e., 

Ve — V = 1 — V , V > Vc (26) 

For V < Vc, we find from geometry Ve = I — f^al + ^(ag, ci|, ■ ■ ■)• Furthermore, the 



amplitude 02 is given by 02 ~ c(f c — f ) 2 , where the coefficient c depends on both co and g. 
Thus, immediately below the transition, we have 

64 2 . 64 2 .X 

- f ~ 1 C Vc + [ C - I) V , V ^ Vc . (27) 

135 435 J ^ ^ c \ J 

This equation implies that all the shifts in the apparent reduced volume for different Co meet 
the curve Vg — v = 1 — v with a slope, (^c^ — 1), larger than —1. Inspection of Fig. p 
suggests that the slope is, in fact, positive for small Co- 

Figures ^ and [1^ show the gravitationally corrected (phase) diagrams with the corrected 
data points for vesicles A and C. The uncertainties of the gravitationally corrected points 
include both the original experimental uncertainties and the numerical uncertainties of the 
gravitational shape-energy calculations. Note that the data points are shifted to the left 
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in f , as expected qualitatively. All experimental points are in the (metastable) prolate 
phase, except the "bad" point which still remains above the upper spinodal, unless the 
"anomalous" cluster is removed, as discussed in the last section. Accepting this somewhat 
ad hoc procedure, one may argue that this last stable point at v = 0.890 lies within error 
bars inside the prolate phase after the gravitational corrections have been performed. We 
note, however, that the actual point of budding at Vb = 0.878, which must in principle be 
beneath the spinodal, would still appear to lie slightly in the unstable region, even after 
gravitational and tipping corrections. We may speculate that this apparent inconsistency is 
due to thermal shifts (neglected so far in our treatment), which could be appreciable near 
the spinodal. 

So far, we have discussed stability of vesicle A only with respect to the upper (pear-mode) 
spinodal. We have also checked that the data points of vesicle A fall above the lower limit 
of stability of the prolates, which is an instability towards the oblate phase. In fact, the first 
point of vesicle A is located (including the gravity correction) at {v = 0.937, Cq = 1-6). The 
location of the prolate-oblate transition at the same volume (and g = 2.2) is known to occur 
at Co = (see the previous section), which is comfortably below our data point. Thus, the 
vesicle-A trajectory does (properly) start in the stable-prolate region. On the other hand, 
this first point still appears somewhat out of line with the remaining three points of the 
trajectory (see Fig. 10), which (after gravity corrections) fit quite well to the simple form 
|l^ , rfiQV = const. We may speculate that thermal shifts play a role here, too, near the 
lower spinodal . 



A few comments are in order concerning gravitational corrections to the vesicle-C data. 
The reduced volume of vesicle C is shifted to the left, as expected. Unfortunately, the 
effective spontaneous curvature Cq is not well determined, due to the large experimental and 
numerical errors. (We remind the reader that vesicle C did bud at this location and, thus, 
exhibited large spinodal fluctuations.) This data point appears to be located appreciably 
below the spinodal line. This could be an artifact created by effects not included in the 
mapping (see the discussion at the end of Sec. 0) and/or it could be due to a low activation 
energy for budding near the sphere. 

The upshot of this exploration of gravitational corrections is that gravity does, indeed, 
have a substantial numerical effect, as might be anticipated from the fact that the dimen- 
sionless parameter g is around 2 for vesicles A and C. However, the qualitative (and generally 
encouraging) conclusions of the gravity-free analysis are not changed. 



VII. CONCLUSION 

Previous experiments (e.g., Ref 0) have compared experiment with theory by, in effect, 
exhibiting a set of control parameters {v, cq, mo, a) which lead to theoretical shapes similar 
to those observed in the laboratory. It is important that this exercise can be successfully 
carried through; but, for various reasons, it constitutes far less than a full test of the theory. 

The first problem is that different variants of the theory, ranging from the SC model 
(a = 0) to the AA (bilayer-couple) model (a = oo), all share the same set of stationary 
shapes, so that observation of a shape which can be suitably parameterized only distin- 
guishes models in which the shape is stable from those in which it is not. One would like to 
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be able to measure all the control parameters for a given vesicle and then to verify that a 
vesicle with those control parameters does, indeed, have the observed shape. The difficulty 
is that, while v and a are measurable, Cq and tuq (which enter the shape problem in the 
combination mo) are not. We have surmounted this problem by concentrating on the equiv- 
alent variable cq (which incorporates cq, mo, and a) and inferring this variable directly from 
the shape data (04). Although the inference process uses theory, there arc nontrivial checks 
left over. Local stability is still an important check, as we have argued. In particular, the 
observed (reasonable) agreement of the experimental budding boundary with the calculated 
theoretical spinodal is encouraging, as is the qualitative agreement of the post-budding shape 
with that predicted by the theory. Once the mapping is done, observation of other shape 
coefficients (02, a^-, etc., for the prolatcs) provides, in principle, a further test of agreement 
between theory and experiment. Unfortunately, at the level of precision we have been able 
to achieve here, 02 is too weakly dependent on cq to be useful, and Og is too small. 

The second and in many ways more important advance which our experiment makes over 
previous ones is in the monitoring and analysis of the full thermal shape ensemble. Previous 
workers have certainly observed the shape fluctuations; however, shape comparisons between 
theory and experiment have heretofore relied on comparison of a single judiciously selected 
image with a theoretically calculated shape. When fluctuations are appreciable (which they 
certainly become near any instability boundary), this process is clearly unacceptable. We 
have illustrated how to monitor and to analyse the full shape ensemble, and we have shown 
how to relate the ensemble data to the corresponding T = theoretical shapes, at least 
in situations where fluctuations are not too large. Treatment of larger fluctuations, which 
are common near instabilities and will certainly be increasingly important at low v (where 
mechanical modes will tend to be softer), will require a new theoretical approach capable of 
going beyond the Gaussian level. 

Finally, at a somewhat technical level, we have illustrated that the effects of gravity, 
which have been ignored in earlier work, are quantitatively important. And, we have shown 
how to adjust for them in comparing theory and experiment. 

In summary, our work provides in principle a quantitative test of the ADE model of 
vesicle shapes. Agreement between theory and experiment (including suitable corrections) 
is crude but satisfactory. It is important at this point (and entirely feasible) to carry out 
similar analyses in other parts of the phase diagram. When more precision becomes available 
in future experiments, more consistency checks will be possible (e.g., by looking at (02) and 
(oe)), and it will be worthwhile to include in the analysis the corrections of order /c^T/o;, 
which we have ignored herein. It is clear that gravitational corrections will have to be 
included and that non-Gaussian effects will be important near instabilities. 
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FIGURES 



FIG. 1. ADE phase diagram for a = 1.4, prolate region. First-order boundaries (D) are 
indicated by solid lines; second-order boundaries (C), by dashed lines; and, spinodals (M), by 
dotted lines. Lowest-energy shapes are illustrated for each region. All symmetry axes are vertical 
as indicated for the prolate. Prolate shapes are locally stable between the upper spinodal line 
Mjj™) and the lower spinodal hue Mg™). The hues DP^/p^"'^, Qpro/pear^ jjpro/obi^ j^pro/nas 

bound the region where prolates are the lowest-energy shapes. In the region immediately above 
-energy shapes are pear-like. LP^"*", is the limiting line at which the neck of the pear 
shape shrinks to zero radius producing a vesiculated shape, as indicated. The lowest-energy states 
above LP^'*'" are dominated by vesiculated shapes. In the region immediately below the prolates, 
oblate and non-axisymmetric (nas) shapes have lowest energy [Q. The point T on the prolate/pear 
phase boundary is a tricritical point, separating first-order and second-order behavior. CEP labels 
a critical end point, where a second-order boundary C^'^^/"^^ (not shown, since it is very close to 
jjpro/nas^ disappears underneath the lower prolate boundary. Note for future reference that the 
limiting line L^'^"''^ crosses the upper spinodal line M^q°_^ of the prolate phase at v = 0.875 for 
a = 1.4. 



FIG. 2. The shape coefficient function of reduced volume v for various values of the 

effective reduced spontaneous curvature cq, as indicated on the curves. These curves allow us to 
infer a value of cq, if a^'^ is known at given v. This is the basis of the mapping procedure discussed 
in Sec. II C| . All the curves pass through a^^^ = at the sphere, v = I. Although 04^^^ is almost 
independent of cq at v = 0.85, there is no common crossing point. 



FIG. 3. Time-sequence of phase-contrast video images of vesicle A at f = 0.954. The images 
are ordered in time from the upper left to the lower right. The elapsed time between images 
is 6.3 s. The length of the long vesicle axis is approximately 20 /im. The vesicle fluctuates 
about an axisymmetric prolate shape; however, each particular contour is different and, in general, 
non-axisymmetric. 



FIG. 4. Time-sequence of phase-contrast video images of vesicle A at v = 0.912. Times and 
scale are as in Fig. |3[ The vesicle is now more elongated than it was in Fig. ^. Strong pear-like 
fluctuations in each direction are now clearly visible. 
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FIG. 5. Time-sequence of phase-contrast video images of vesicle A {Ra = 9.2 /um) at u ~ 0.878, 
illustrating the budding process. The scale is as in Fig. ^ and |^. The time elapsed between images 
here is 1.2 s. A pear-like fluctuation, much like those visible in Fig. |^, now carries the vesicle over 
the metastable barrier to the budded state. The pear shaped contours correspond to transient 
shapes and are not stable. Note that the ratio of the vesicle size to the satellite size after budding 
is roughly 2.8. 



FIG. 6. Time-sequence of phase-contrast video images of vesicle B {Ra = 5.5 fj,m) at u ~ 0.945, 
illustrating the budding process for a smaller size vesicle. The time elapsed between images is the 
same as in Fig. ^. Here, the transition from the prolate via the transient pear to the budded state 
happens much more quickly than for vesicle A, due to the smaller hydrodynamic radius. 



FIG. 7. Typical time-series for the amplitudes as and 04 and for the effective reduced volume Vg 
for vesicle A at reduced volume (ve) = 0.954. The dashed line corresponds to the mean amplitude 
(an), which is close to zero for as ||5^. According to Eq. (|2^), the mean values {ve) and (a4) 



correspond to the zero-temperature quantities v and 04^^^ , respectively, which are the basis for the 
mapping. 



FIG. 8. Same data as in Fig. |^ but for reduced volume (ve) = 0.912. Note the longer 
time-scales for fluctuations in the as-mode as the shape instability (spinodal) is approached. 
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FIG. 9. Experimental trajectories without (Sec. V, open symbols) and with (Sec. VI, filled 
symbols) gravitational corrections in the {v,co) diagram for vesicles A, B, and C. Corrected and 
uncorrected data points are joined by a thin line. As explained in the text, the vesicle follows 
a path from lower right to upper left, as it is heated. Vesicles A and B underwent a budding 
instability after the upper-left-most point of the trajectory, vesicle C budded from the point shown. 
The spinodal lines MP^° are the same as those shown in Fig. || and mark the upper and lower 
boundaries of the region predicted by theory to be locally stable for prolate shapes. Thus, for full 
consistency, the trajectories must terminate below the upper spinodal. Vesicles B and C satisfy 
this criterion. The final raw-data point on the trajectory of vesicle A is inconsistent with stability. 
Including gravitational eff'ects moves the thermal trajectories to smaller reduced volume. For 
g = 2.2 (appropriate only for vesicle A), the instability boundaries are shifted by gravity, as shown. 
The final point on the vesicle A trajectory becomes consistent with theory, only when, in addition 
to gravitational effects on the mapping, the gravitational tipping of the fiuctuating pear-like shapes 
is incorporated by excluding a cluster of data points (square symbols) , as discussed at the end of 
Sec. 0. Uncertainties are generally large near the spinodal lines and close to the sphere, where 
fluctuations become important. 



FIG. 10. Experimental trajectories without (Sec. V, open symbols) and with (Sec. VI, 
filled symbols) gravitational corrections in the (vjItiq) phase diagram for vesicles A, B, and C 
using a = 1.4 (see Fig. |9| for a legend). The stability of data points does not depend on a 
and mirrors Fig. |^. As for Fig. |^, the last point of the vesicle A trajectory becomes consistent 
with the theoretically calculated stability when corrected for gravitational effects (including the 
cluster exclusion, as discussed in the text). A theoretical thermal trajectory (with the simple 
form fhov = const) is shown for comparison. Although Figs. ^ and [l^ look similar, they are not 
connected by a simple rescaling of the vertical axis. 



FIG. 11. Computed gravitational corrections to the effective reduced volume Ve at g = 2.2, 
as is appropriate for vesicle A. The shift {ve — v) is shown vertically as a function of cq and v. 
The dashed line corresponds to the prolate-oblate transition. All shifts meet this line with a slope 
different from —1, as explained in the text. The shaded region (fg > 1) is unphysical. 
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TABLES 



Vesicle 


Tj C 


/„, \ 


(04) 




Z6.7 ± U.Ui 


u.yo4 ± u.uui 


A AAO'7 I A AA1 1 

— 0.0U37 ± U.UOll 




0Z.7 ± U.UI 


u.yoz ± U.UUI 


U.UU54 ± U.UU13 


A 


37.8 ± U.UI 


U.912 ± U.OUl 


0.0U38 ± 0.0U22 




42.4 ± U.UI 


0.894 ± O.OUl 


U.U092 ± 0.001 




45.8 ± U.UI budding 


U.878 ± 0.002 


N/A 




32.9 ± U.UI 


0.970 ± 0.001 


0.005 ± 0.001 


B 


4U.3 ± U.UI 


0.950 ± O.OUl 


U.U16 ± O.OUl 




42.U±U.U1 budding 


0.945 ± 0.0U2 


N/A 


C 


27.0 ± 0.01 budding 


0.983 ± 0.001 


0.008 ±0.002 



TABLE I. Experimental results for vesicles B, and C. At each temperature T, the values 
of {vg), and (04) are shown. The final budding temperature is also given, along with the (for A 
and B extrapolated) reduced volume at budding. 
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Fig. 9: Dobereiner et al. 
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Fig. 10: Dobereiner et al. 
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